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ABSTRACT 

Close-in exoplanets with highly eccentric orbits are subject to large variations in incoming stellar flux 
between periapse and apoapse. These variations may lead to large swings in atmospheric temperature, 
which in turn may cause changes in the chemistry of the atmosphere from higher CO abundances at 
periapse to higher CH4 abundances at apoapse. Here we examine chemical timescales for CO^CH4 
interconversion compared to orbital timescales and vertical mixing timescales for the highly eccentric 
exoplanets HAT-P-2b and CoRoT-lOb. As exoplanet atmospheres cool, the chemical timescales for 
CO«=^CH4 tend to exceed orbital and/or vertical mixing timescales, leading to quenching. The rela- 
tive roles of orbit-induced thermal quenching and vertical quenching depend upon mixing timescales 
relative to orbital timescales. For both HAT-P-2b and CoRoT-lOb, vertical quenching will determine 
disequilibrium CO«=^CH4 chemistry at faster vertical mixing rates (K zz > 10 7 cm 2 s" 1 ), whereas 
orbit-induced thermal quenching may play a significant role at slower mixing rates (K zz < 10 7 cm 2 
s _1 ). The general abundance and chemical timescale results - calculated as a function of pressure, 
temperature, and metallicity - can be applied for different atmospheric profiles in order to estimate 
the quench level and disequilibrium abundances of CO and CH4 on hydrogen-dominated exoplanets. 
Observations of CO and CH4 on highly eccentric exoplanets may yield important clues to the chemical 
and dynamical properties of their atmospheres. 

Subject headings: planetary systems — planets and satellites: atmospheres — planets and satellites: 
composition — planets and satellites: individual (HAT-P-2b, CoRoT-lOb) 



1. INTRODUCTION 

For planets with high eccentricity, the large vari- 
ations in flux received from their host stars may 
yield substantial variations in atmospheric tempera- 
ture a nd dynamical behav i or du r ing the course of an 
orbit (iLangton fe Laughlm] 120081: lLaughlin et ail 120091: 
Iro fe peming||2010l;lCowan fe Agolll201ll;lKane fe Gelinol 
20111: iRauscher fe Menoul 120121: [Lewis et al.l 12012ft . In 
some cases, orbit-induced temperature variations may 
be large enough to produce a significant shifts in the 
chemical behavior of the planet. For example, the swing 
between high atmospheric temperatures at periapse and 
lower atmospheric temperatures at apoapse may shift 
equilibrium chemistry predictions from relatively higher 
CO abundances at periapse (toward a CO-dominatcd 
atmosphere) to relatively higher CH4 abundances at 
apoapse (toward a CH4-dominatcd atmosphere). 

Large changes in the abundances of CO and CH4 are 
of particular interest because these compounds strongly 
influence the spectral properties of e x oplanet atmo- 
spher e s (e.g.. ISeager fe Sasselovl 120001: IBurrows et al 
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Lewis et al.l |2012|) may 
in chemical composition of highly eccentric transiting 



exoplanets (see Table Q]), part icularly at wave lengths 
sensitive to CO and CH4 (e.g. JLewis et al.ll2012ft 

However, the extent of temperature-dependent varia- 
tions in carbon chemistry throughout the orbit - and 
whether equilibrium chemistry at a given altitude pre- 
vails over orbital timescales - also depends upon the rate 
of CO«=^CH4 interconversion relative to the time elapsed 
between periapse and apoapse. Thermochemical equi- 
librium can be maintained throughout the orbit only 
if chemical timescal es are less than orb ital timescales. 
To study this effect. Uro fc Demingl (|2010f ) estimated the 
CO^CH4 interconversion time scale for HP 80606b and 
HD 17156b using the kinetics of Bezard et al. (2002]) and 
found that orbital timescales are generally much shorter 
than chemical timescales - indicative of disequilibrium 
chemistry - at pressure levels where orbit-induced tem- 
perature variations are expected to be significant. This 
behavior will occur on objects which are expected to 
have relatively low atmospheric temperatures (and there- 
fore sluggish reaction kinetics) near apoapse, leading to 
orbit-induced thermal quenching wherein an equilibrium 
composition achieved near periapse survives to become a 
disequilibrium composition at apoapse. 

The observed properties of highly eccentric exoplan- 
ets are subject to numerous variables including ther- 
mochemical and photochemical reaction rates, convec- 
tive transport, and horizontal dynamical and radiative 
timescales. Here, we focus specifically on the role of 
thermochemical quench chemistry in response to vertical 
transport and eccentricity-induced atmospheric temper- 
ature variations, as this behavior may strongly influence 
observable CO and CH4 abundances (even if photoly- 
sis occurs at higher altitudes). For simplicity, the rele- 
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TABLE 1 

Highly-Eccentric (e > 0.3) Transiting Exoplanets 



Object 


o(AU) 


e 


p(days) 


Ps (days) 


M 3 


Rj 


HD 80606b 


0.447 


0.934 


111.44 


1.7 


3.9 


1.03 


HD 17156b 


0.163 


0.682 


21.22 


3.6 


3.3 


1.02 


CoRoT-lOb 


0.105 


0.530 


13.24 


4.3 


2.8 


0.97 


HAT-P-2b 


0.068 


0.517 


5.63 


1.9 


8.9 


1.16 


HAT-P-34b 


0.068 


0.440 


5.45 


2.4 


3,1 


1.20 


HAT-P-17b 


0.088 


0.346 


10.34 


5.9 


0.5 


1.01 


WASP-8b 


0.080 


0.310 


8.16 


5.1 


2.1 


1.04 


Note. 


Data from 1 Wright et al.l ( 


2TJlT|). 


The pseudo- 



synchro nous rotati on period p„ wa s calculated using; the expres- 
sions of lHull (QUID as presented in llro &; De ming (201^). 

vant timescale for the temperature variation is taken to 
be 0.5p (where p is the orbital period), which describes 
the time elapsed between temperature swings at periapse 
and apoapse. We first calculate the abundances of CO 
and CH4 and CO^CHi chemical timescales as a func- 
tion of pressure, temperature, and metallicity in a solar- 
composition ga s, using updated kinetics for CO<=^CH4 
interconversion (IVisscher et al.ll2010bl:lMoses et al.N2011fc 
IVisscher fc Mosesl 1201 ID ~ The results are then com- 
pared to orbital and mixing timescales and pressure- 
temperature profiles of individual objects to estimate the 
quench levels and abundances of CO and the CH4 as dis- 
equilibrium species. The primary objective of this study 
is to discuss a convenient method for exploring ther- 
mal quenching processes in the CO^CILi system in exo- 
planet atmospheres. Although we focus on CO-abundant 
HAT-P-2b and CH 4 -abundant CoRoT-lOb as specific ex- 
amples, the abundance and timescale results presented 
here may in principle be applied to any H2-dominated 
substellar object that is subject to orbit-induced tem- 
perature variations. 

2. CHEMICAL MODEL 

2.1. Equilibrium Chemistry and Abundances 

We develop a gas-phase chemical mo del based upon the 
therm oc hemical equilibrium mode ls of iFeglev fc Lodd"eri 



199 1 lLodders fc Feglevl (I2002T). and IVisscher et al 
200a I2010af) . As in lLodders fe Feglevl (|2002|), we ex- 
amine H-C-N-0 chemistry over a range of pressures, 
temperatures, and metallicities so that the model re- 
sults may be applied to a variety of substellar objects. 
Solar-composit ion elemental abundances are taken from 
lLodders! ((2010), and we consider the re moval of ~ 21% 
of oxygen into rock-fo rming elements (lLodders ! 120041 : 
IVisscher fe Feg lcv 2005). Figure [T] shows the equilibrium 
mole fraction abundances of CO and CH4 as a function 
of pressure and temperature in a solar-metallicity gas. 
Also shown is the CO=CH4 equal abundance boundary 
(Ago = A C h 4 = 0.5A SC ), approximated by 



logP w 5.05 - 5807.5/T + 0.5[Fe/H], 



(1) 



for P in bars and T in kelvins. At pressures and tem- 
peratures where methane is the dominant C-bearing gas 
(higher P, lower T), the mole fraction abundances of 
CH4 and CO are given by 



Ach 4 
Xco 



A£ C 771, 

A co (P,T)m 2 



(2) 
(3) 



where X^ c is the total carbon abundance in a solar- 
metallicity gas (~ 4.5 x 10- 4 ; ILodderl[20T0t) . A co (P,T) 
is the solar-metallicity abundance of CO as a function of 
pressure and temperature as plotted in Figure [T^, and 
i ji is the metallicity factor defined by log m = [Fc/H] 
(jLodders fc Feglevl 12001 IVisscher et atl l2006). 

At pressures and temperatures where carbon monoxide 
is the dominant C-bearing gas (lower P, higher T), the 
mole fraction abundances of CO and CH4 are given by 



Ago 



X£ c m, 



(4) 
(5) 



where A CH4 (P, T) is the mole fraction abundance of CH4 
as a function of pressure and temperature in a solar 
metallicity g plotted in Figure [Tp. Unlike CO, 

which shows a strong (m 2 ) dependence on metallicity in- 
side the CH4 field, the methane abundance inside the CO 
field is independent of metallicity because it is propor- 
tional to the total carbon abundance (SC) but inversely 
proportional to the total oxygen abundance (SO). Note 
that these simple abundance approximations begin to 
break down at higher metallicities (m > 30) as heavy el- 
ements become significantly abundant relative to hydro- 
gen and helium. For high-metallicity cases, mole fraction 
values can be calcula ted from solar elemental abundance 
ratios ()Lod dcrs 2010) with different values of m. A more 
complete description of the equilibrium chemical behav- 
ior of CO, CH4, and other carbon-bearing species as a 
function of P, T, and m in substell ar atmospheres can 
be found in lLodders fc Feglevl (|2002fl . 

2.2. Chemical Kinetics and Timescales 

For CO- CH4 quench k inetic s we a dopt the reaction 
schem e of iMoses et al.l (|2011| ) and IVisscher fc Mosesl 
(|2011[) . who identify two plausible mechanisms for 
CO«=^CH4 interconversion in substellar atmospheres. In 
general, the rate-limiting step for breaking/forming the 
C-0 bond is 

CH3OH + M CH 3 + OH + M, (6) 

when CH4 is the dominant C-bearing gas, and 

CH2OH + H ^ CH 3 + OH, (7) 

when CO is the dominant C-bearing gas. However, 
the contribution of either pathway is significant enough 
that both should be considered when estimating the 
CO^CH4 interconversion rate in substellar atmospheres. 
The chemical timescale for CO— >CH4 is thus given by 



Tchem(CO) 



[CO] 



^[CH 3 OH][M] + / Q Cll,()ll II ' 



(8) 



where Ajgj^- and ^yj^ are the forward reaction rate con- 
stants for reactions (|6]) and ([7]), respectively. The chem- 
ical timescale for CH4— >CO is given by 



7"chem(CH 4 ) — 



[CH 4 



^[CH 3 ][OH][M] + ^[CH 3 ][OH] ' 



(9) 



where and ^jjj, are the reverse reaction rate con- 
stants for reactions ([6|) and (0, respectively, taken from 
IJasper et al.l ()2007D and reve rsed to give and Iff^ 

using the method described in IVisscher fc Mos es (2011). 
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Fig. 1. — Equilibrium mole fraction abundance contours (on a logarithmic scale) for (a) CO and (b) CH4 in a solar-composition gas 
at pressures and temperatures relevant for extrasolar giant planet atmo spheres. In each panel, the dashed line denotes the CO = CH4 
equal-abundance boundary where Xco = XrjH 4 ~ 0-5X^c- Figure after Loddcrs & Fcglcv (2002) 



Using equilibrium abundances for OH, CO, CH 3 , 
CH4, CH2OH and CH3OH, the chemical timescales for 
CH4— ^CO and CO— ^CHi are shown in Figure [5] (as con- 
tours on a logarithmic scale) as a function of pressure 
and temperature in a solar-metallicity gas. We also 
calculated chemical timescales at metallicities between 
[Fe/H] = -2 and [Fe/H] = +2. Over this range, the 
timescale for CO destruction, r c h e m(CO), is independent 
of metallicity, whereas the timescale for CH4 destruction 
as a function of metallicity is given by 

Tchcm(CH 4 ) « r* hom (CH 4 )(P, T)m-\ (10) 

where 7" c *h om (CH4)(P, T) is the chemical timescale for 
CH4— >CO as a function of pressure and temperature 
in a solar metallicity gas (plotted in Figure [5Jd) and in 
is the metallicity factor. Equation (|10p demonstrates 
that higher metallicities result in shorter timescales (i.e., 
faster reaction kinetics) for CH4— >CO chemistry over the 
range of pressures and temperatures considered here. 

3. RESULTS AND DISCUSSION 

The abundances in Figure [1] and the chemical 
timescales in Figure [5] can be compared with orbital and 
dynamical timescales as well as atmospheric pressure- 
temperature profiles for individual objects in order to 
estimate the quench level and abundances of disequilib- 
rium species. Using CO-dominated HAT-P-2b and CH4- 
dominated CoRoT-lOb as specific examples, here we ex- 
amine various profile/timescale intersections in order to 
explore C0<=±CH4 quench processes in exoplanet atmo- 
spheres, as these intersections represent some quench 
level with respect to orbital timescales. The guiding 
principle throughout the following discussion is that the 
chemical behavior in the CO^CHi system is governed 
by the shortest available timescale in the system. 

The planet-wide average temperature profiles at pe- 
riapse and apoapse overlaid on the timescale and abun- 
dance plots are shown in Figure [3] for HAT-P-2b and Fig- 



ure |4] for CoRoT-lOb. The pressure-temperature profiles 
for HAT-P-2b are planetary average profiles derived from 
the general circulation models (GCMs) of iLewis et al.1 
(|2010Ll2012fl . The Co RoT-lOb profile s are f r om the one- 
dimensional models of iFortnev et al.l (120061 120081 12010ft 
where the incident stellar radiation is assumed to be re- 
distributed over the entire planet. 

Because we are primarily interested in the temperature 
difference between periapse and apoapse (and its effect 
on thermochemistry), we focus on the time elapsed be- 
tween these extremes irrespective of radiative timescales. 
However, orbit-induced temperature variations will ex- 
ist only if the radiati ve timescale is less th an the or- 
bital period (e.g., see iCowan fc Agol 1201 ll and refer- 
ences therein). Orbital timescales for the highly eccen- 
tric exoplanets listed in Table Q] are on the order of days, 
whereas radiative timescales at photospheric altitudes on 
such planets a r e typ i cally on the order o f hours (e.g., 
Laughlin et all 120091: llro fc Demine] 120101: ILewis et al.l 
2012h . We thus adopt 0.5p (where p is the orbital period) 
as the characteristic orbital timescale throughout the fol- 
lowing discussion, and time elapsed between periapse and 
apoapse is denoted in Figures [3] and 0] by the lines la- 
beled 0.5p. The time elapsed between noon and midnight 
is denoted by 0.5p s , which is one half ofthe pseudo- 
synchronous rotation period deriv ed from iHutJ (|1981l ) as 
presented in llro fc Demind (|2010h . Although we focus 
on orbital temperature variations in the present study, 
in principle diurnal variations in atmospheric chemistry 
could be explored using a similar approach, depending 
upon the relative dynamical and radiative timescales. 



3.1. Transport-induced Quenching 

In addition to orbit-induced temperature variations, 
atmospheric constituents will also be subject to convec- 
tive vertical mixing. The mixing timescale (r m i x ) is esti- 
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Fig. 2. — Chemical timescales (in seconds, on a logarithmic scale) for (a) the conversion of CO to CH4 and (b) the conversion of 
CH4 to CO in a solar-m etallicity gas. The dashed line denotes the position of the CO = CH4 equal abundance boundary. Figure after 
ILodders fc~ Fcglcv (2002?). See text for details. 



mated using the expression 

r mix = L 2 /K zz , (11) 

where L is the mixing length and K zz is the eddy diffu- 
sion coefficient (|Stonelll976t ISmithl Il998). Because the 
vertical mixing rate in extrasolar giant planets is un- 
known, K zz is treated as a free parameter in the models. 
For hot giant plan et atmospheres, it is rea s onable to as- 
sume L ~ 0.5H (jVisscher fc Mosesl 120111: iSmithl 11998ft 
where H is the atmospheric scale height (H = RT/fig). 
If required, better estimates of L/H may b e calcu l ated o n 
a case-by-case basis using the method of ISmithl ([1998). 
The gravity is calculated from the mass and radius of 
the planet (see Table [IJ and we assume a mean molec- 
ular mass of fx ~ 2.4 g mole -1 for an ^-dominated, 
solar-composition atmosphere. 

The vertical quench level i s defined as the alti- 
tude at which T m j x = T c h om l|Prinn fe B arshavj 119771 ; 
iFeriev fc Pri^l[l985h . The pressures and temperatures 
for which this condition occurs on each planet are shown 
in Figures [3] and |4] as filled circles with labeled solid 
lines denoting the estimated K zz value (in cm 2 s -1 ) for 
each case. The intersections of the atmospheric pressure- 
temperature profiles with the r m ; x = r c h e m curves in Fig- 
ures [3] and [4] indicate the location of transport-induced 
quench level (P q , T q ) for different assumptions about the 
vertical mixing rate (characterized by K zz ). The quench 
level represents the lowest pressure (i.e., highest altitude) 
along the atmospheric profile for which equilibrium can 
be achieved in the presence of convective vertical mixing; 
at all pressures < P q , the abundance of CO (for CO de- 
struction) or CH 4 (for CH 4 destruction) will remain fixed 
at the equilibrium abundance achieved at P qi T q . As ex- 
pected, high K zz values result in quenching deep in the 
atmosphere whereas low K zz values give quenching at 
lower atmospheric pressures. However, as discussed be- 
low, vertical mixing can only play a significant role in the 
atmospheres of highly eccentric exoplanets over orbital 



timescales if mixing timescales (T m ; x ) are shorter than 
t he time elapsed betwe e n per iapse and apoapse (0.5p). 

ICooper fc Showman! (|2006| ) demonstrated that dise- 
quilibrium abundances of CO and CH4 in the upper (ob- 
servable) atmospheres of hot Jupitcrs arc controlled by 
vertical quenching, and that horizontal transport may 
homogenize these abundances with longitude. In the 
present study we therefore do not examine the effects 
of horizontal transport, and focus instead on the rela- 
tive roles of vertical quenching caused by convective mix- 
ing and thermal quenching caused by orbit-induced tem- 
perature variations. However, we note that strong hor- 
izontal winds may serve to homogenize disequilibrium 
abundances from either quenching mechanism (vertical 
or orbit-induced) in the upper atmospheres of highly ec- 
centric exoplanets. 

3.2. Methane Quenching on HAT-P-2b 

On HAT-P-2b, CO is the dominant carbon-bearing gas 
and so we consider CH4 quench chemistry. As the planet 
cools between periapse and apoapse, the equilibrium CH4 
abundance is predicted to greatly increase (e.g., see Fig- 
ure [3];). However, equilibrium can only be maintained if 
CO— >CH4 conversion kinetics are fast enough to produce 
CH4 over orbital timescales (r c h C m < 0.5p). We therefore 
examine CO— J>CH4 chemical timescales (left hand panels 
in Figure |3|) in order to determine whether the higher 
CH4 abundances predicted by equilibrium can be pro- 
duced as the atmosphere shifts from high temperatures 
at periapse to low temperatures at apoapse. 

At relatively high pressures and temperatures, deep in 
the atmosphere, equilibrium can be maintained through- 
out the entire orbit if chemical timescales are always 
shorter than orbital timescales (r c h C m < 0.5p). For 
highly eccentric exoplanets, the lowest pressure for which 
this condition holds is given by the intersection of the 
apoapse profile and the orbital timescale (0.5p). On 
HAT-P-2b, the apoapse/0.5p intersection occurs near the 




Fig. 3. — Chemical timescales (top panels; contours in seconds on a logarithmic scale) and equilibrium CH4 abundances (bottom panels; 
mole fraction contours on a logarithmic scale) for CO— >CH4 conversion timescales (left panels) and CH4— >CO conversion timescales (right 
panels) on HAT-P-2b. Planet- wide average temperature profiles for HAT-P-2b are shown for apoapse (blue lines) and periapse (red lines). 
The solid green lines denote the time elapsed between periapse and apoapse (0.5p, where p is the orbital period) and the time elapsed 
between noon and midnight (0.5p s , where p s is the pseudo-synchronous rotation rate). The filled symbols with tie lines show where 
T mix = T cho m the vertical quenc h level) for different values of K zz (cm 2 s — 1 ) using the mixing length approach of Smith (1998) as 

presented in Visschcr & Moses (2011). (A color version of this figure is available in the online journal.) 



4-bar level (see Figured), demonstrating that CO— >CH 4 
equilibrium chemistry - and the corresponding increase 
in the equilibrium abundance of CH 4 - can be main- 
tained throughout the entire orbit only for P > 4 bars. 
In the absence of vertical mixing, the equilibrium CH 4 
mole fraction abundance at the ~ 4-bar level can track 
eccentricity-induced temperature variations over orbital 
timescales by shifting from ~ 1 ppm at periapse to ~ 10 
ppm near apoapse (see Figure [3};). 

At lower pressures, however, orbit-induced tempera- 
ture variations will quench CO— >CH4 conversion once 
Tchem > 0.5p as the planet moves between periapse 



and apoapse. Equilibrium compositions achieved dur- 
ing part of the orbit when the atmosphere is warmer 
(T c hem < 0.5p) will quench once Tchem > 0.5p as the at- 
mosphere cools. This boundary is represented by the 0.5p 
curve between the apoapse and periapse profiles in Fig- 
ure Quenched CH4 abundances resulting from orbit- 
induced thermal quenching can likewise be estimated by 
the position of the orbital timescale (0.5p) curve in Fig- 
ure Note that the portion of the orbit over which 
equilibrium is achieved (i.e., T c hem < 0.5p) decreases with 
altitude in the upper atmosphere of HAT-P-2b. 
The intersection between the periapse profile and the 




Fig. 4. — Chemical timescales (top panels; contours in seconds on a logarithmic scale) and equilibrium CO abundances (bottom panels; 
mole fraction contours on a logarithmic scale) for CO— >CH4 conversion timescales (left panels) and CH4— >CO conversion timescales (right 
panels) on CoRoT-lOb. Planet-wide average temperature profiles for CoRoT-lOb are shown for apoapse (blue lines) and periapse (red 
lines). The solid green lines denote the time elapsed between periapse and apoapse (0.5p) and the time elapsed between noon and midnight 
(0.5p s ). The filled symbols with tie lines show where T m ; x = T c h cm (i.e., the vertical quench level) for different values of K zz (cm 2 s -1 ). 
(A color version of this figure is available in the online journal.) 



0.5p curve represents the highest altitude (i.e., lowest P) 
in the atmosphere at which equilibrium can be achieved 
at any point during the orbit, with respect to orbital 
timescales. In Figure |3ji, the intersection between pe- 
riapse and the 0.5p curve occurs near the 1-bar level. 
At this intersection, CO— >CH 4 proceeds over orbital 
timescales because r c h cm (CO) = O.bp. As temperatures 
cool from periapse to apoapse (moving left horizontally 
from the periapse/0.5p intersection in Figure [3^), the 
rate of CO destruction (and hence CH 4 production) de- 
creases and quenches when r c hem(CO) > 0.5p. For HAT- 
P-2B, this suggests orbit-induced thermal quenching of 
~ 0.1 ppm CH4 at the ~ 1-bar level throughout the or- 



bit, far below the equilibrium abundance of ~ 10 ppm 
predicted at apoapse temperatures (see Figure Ht). 

However, this thermal quenching scenario neglects 
the role of vertical transport. The relative importance 
of vertical versus orbit-induced thermal quenching de- 
pends upon the transport timescale relative to the or- 
bital timescale. As can be seen in Figure [3J the or- 
bital timescale for HAT-P-2b is roughly equivalent to 
a mixing timescales characterized by K zz ~ 10 7 cm 2 
s _1 . For K zz < 10 7 cm 2 s _1 , orbit-induced thermal 
quenching from periapse to apoapse will likely determine 
the CH4 abundance because vertical mixing timescales 
exceed the time elapsed between periapse and apoapse 



Chemical Timescales on Highly Eccentric Exoplanets 
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(jmix > 0.5p). Thermal quenching in this scenario is ex- 
pected to yield very low (0.01 — 1 ppm) CH4 abundances 
in the upper atmosphere throughout the orbit. 

If we assume faster vertical mixing rates, vertical 
quenching will dominate because T c h em < 0.5p. For 
K zz > 10 7 cm 2 s _1 , transport-induced quenching will 
give ~ 10 ppm CH4 throughout the upper atmosphere of 
HAT-P-2b at apoapse (see Figure^) , because thermally- 
quenched CH4 will be overwhelmed by much larger quan- 
tities of CH4 transported from higher pressures alo ng the 
apoapse profile. As noted bv iLewis et al.l (|2012f ). even 
higher CH4 abundances may be possible if HAT-P-2b 
has a red uced metallicity an d/or an enhanced C/O ratio 
(e.g., see lMoses et al.ll2012D . Estimates of the quenched 
CH4 abundance and the expected quench mechanism are 
summarized in Tabled The large difference in predicted 
abundances as a function of K zz suggest that observa- 
tions of the CH4 abundance on highly eccentric exo- 
planets may yield information about atmospheric mixing 
rates. 

3.3. Carbon Monoxide Quenching on CoRoT-lOb 

On the cooler CoRoT-lOb, CH4 is the dominant C- 
bearing gas and so we consider CO quench chemistry. 
Unlike CH4 chemistry on HAT-P-2b, the equilibrium CO 
abundance on C0R0T-IOB is expected to decrease as the 
planet moves from periapse to apoapse (see Figure 
leading to quenching of excess CO. We thus first con- 
sider CO— >Cti4 kinetic timescales as the planet nears 
apoapse. The intersection of the apoapse profile and or- 
bital timescale (0-5p) occurs near the 90-bar level in Fig- 
ure I3J:, indicating that CO— ^CELi equilibrium chemistry 
can be maintained throughout the orbit if P > 90 bars. 

As the planet swings back toward periapse, the equi- 
librium abundance of CO is predicted to rapidly increase 
with increasing temperatures. However, CO can only 
be produced if CH4— >CO kinetics (Figure Hb) are suf- 
ficiently fast relative to orbital timescales. We must 
therefore examine CH4— >CO chemical timescales to de- 
termine the highest altitude (lowest P) at which equilib- 
rium can be maintained throughout the orbit. The inter- 
section between the apoapse profile and the 0.5p curve 
for CH4— >CO timescales occurs near the 300-bar level on 
CoRoT-lOb: i.e., only for atmospheric pressures greater 
than ~ 300 bars can CO^CH^ equilibrium chemistry be 
maintained over orbital timescales, shifting from ~ 0.3 
ppm CO at apoapse to ~ 3 ppm CO at periapse (see 
Figure HJi). However, note that CO abundances in the 
upper atmosphere of CoRoT-lOb at apoapse will more 
likely be determined by orbit-induced thermal quench- 
ing or vertical quenching, depending upon the rate of 
vertical mixing in its atmosphere. 

The intersection between periapse and the 0.5p curve 
represents the lowest pressure level at which equilibrium 
can be achieved at any point during the orbit with re- 
spect to orbital timescales. For CoRoT-lOb, the inter- 
section between the periapse profile and the 0.5p curve 
occurs near 14 bars for CO— !>CH4 chemistry (Figure 
At pressure levels between the apoapse/0.5p inter- 
section (~ 90 bars) and the periapse/0.5p intersection 
(~ 14 bars), CO— >CHi conversion is quenched whenever 
Tchem > 0.5p (this boundary is represented by the 0.5p 
curve between the apoapse and periapse profiles in Fig- 
ure^). From a comparison of orbital timescales (0.5p 



TABLE 2 

Estimated Quench Abundances at P < 1 bar 



K zz 


HAT-P-2b 


CoRoT-lOb 


(cm 2 s" 1 ) 


CH 4 


CO 


10 3 


0.01-0.1 ppm {t) 


>100 ppm {t) 


10 s 


0.03-0.1 ppm (t) 


>100 ppm (t) 


10 7 


10 ppm (v) 


100 ppm (v, t) 


10 9 


10 ppm (v) 


10 ppm (v, t) 



Note. — The notations (v) and (t) denote 
whether the quenching mechanism is predomi- 
nantly caused by vertical mixing or by orbit- 
induced thermal variations, respectively. 



curve) and CO abundances in Figure |4j:, we thus might 
expect thermally-quenched CO abundances of ~ 0.1 — 100 
ppm above the ~ 90-bar pressure level on CoRoT-lOb. 

For example, orbit-induced thermal quenching at the 
~ 14-bar level (the periapse/0.5p intersection) would 
yield ~ 100 ppm CO throughout the orbit, in excess of 
the equilibrium abundance of ~ 1 ppm CO predicted 
at apoapse. This will also overwhelm the ~ 0.1 — 1 
ppm CO predicted from vertical quenching over a range 
of K zz values (10 3 — 10 7 cm 2 s _1 ). Once present at 
lower pressure levels, CO cannot readily be converted 
back into CH4 because r c hom(CO) > 0.5p (i.e., it is 
quenched) and r c h em (CO) > r c h G m(CH4) (i.e., the CO 
production rate exceeds the CO destruction rate). Al- 
though CO production via CH4— >CO is also quenched 
at these altitudes, high abundances of CO are main- 
tained because there is no way to destroy it thermo- 
chemically. In addition, as noted above, rapid horizontal 
transport would tend to homogenize this d isequilibrium 
abundance (e.g.. iCooper fc Showman 20Ql). For slower 
mixing rates (K zz < 10 7 cm 2 s _1 ) it thus appears plau- 
sible that thermal quenching from eccentricity-induced 
temperature variations determines the CO abundance in 
the upper atmosphere of CoRoT-lOb. 

As for HAT-P-2b, the orbital timescale for C0R0T- 
10b is roughly equal to vertical mixing timescales char- 
acterized by K zz ~ 10 7 cm 2 s _1 , and we therefore expect 
vertical mixing to play a role for K zz > 10 7 cm 2 s _1 . 
Note that for these cases, vertical quenching along the 
periapse profile may determine a CO abundance which 
subsequently undergoes thermal quenching at lower pres- 
sures. For example, quenching at K zz ~ 10 7 or ~ 10 9 
cm 2 s^ 1 at periapse would yield quenched CO abun- 
dances of ~ 100 or ~ 10 ppm, respectively, throughout 
much of the upper atmosphere. These quantities of CO 
would prevail via orbit-induced thermal quenching as the 
planet moves from periapse to apoapse and overwhelm 
disequilibrium CO abundances that result from vertical 
quenching (r c h C m = Tmix) along the apoapse profile. 

Estimates of the upper atmospheric (P < 1 bar) 
quenched CO abundance and the expected quench mech- 
anism^) are summarized in Table O Although the re- 
sults for CoRoT-lOb do not show large differences in 
abundance as a function as K ZZl they consistently pre- 
dict disequilibrium CO abundances 2-4 orders of mag- 
nitude in excess of equilibrium abundances predicted at 
apoapse. 

3.4. Nitrogen Quench Chemistry 
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Visscher 



Other chemical systems may also respond to 
eccentricity-induced temperature variations. For exam- 
ple, large changes in atmospheric temperature may shift 
equilibrium chemistry predictions from relatively higher 
N 2 abundances at periapse to relatively higher NH 3 
abundances at apoapse. To explore this effect, we cal- 
culated chemical timescales for N2«=^NH3 interconver- 
si on using the updat ed reaction mechanism described 
in iMoses et al.l ([201 ID . For highly eccentric exoplanets, 
the quench process of most interest is N2— >NtL3 as the 
planet swings from periapse to apoapose (where higher 
NH3 abundances are predicted). However, compared to 
the CO^CELi system, ^^NELj interconversion gener- 
ally quenches at higher pressures, deeper in the atmo- 
sphere, where orbit-induced temperature variations are 
smaller. We thus expect orbit-induced thermal quench- 
ing to play a relatively minor role on HAT-P-2b and 
CoRoT-lOb, yielding variations in the NH3 abundance 
of less than an order of magnitude. Vertical quenching 
will more likely determine the abundance of NH3 in the 
upper atmospheres of these highly eccentric exoplanets. 



4. SUMMARY AND CONCLUSIONS 

Equilibrium abundances for CO and CH4 and chemical 
timescales for CO^CtLj interconversion were calculated 
as a function of pressure, temperature, and metallicity. 
A comparison of the abundance and timescale plots with 
atmospheric pressure-temperature profiles can be used 
to estimate the quench levels and abundances of CO and 
CH4 for a variety of substellar objects. In principle, this 
approach can be applied for any substellar object with an 
H2-dominated atmosphere that is subject to eccentricity- 
induced temperature variations. Overall, orbit-induced 
thermal quenching tends to favor CO over CH4 because 
CO is the higher-temperature species and because chem- 
ical reactions proceed more rapidly at periapse than at 
apoapse. Moreover, equilibrium can be maintained to 
higher altitudes on warmer CO-dominated objects than 
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on cooler CELi-dominated objects, with respect to orbital 
and vertical mixing timescales. 

Whether vertical quenching or orbit-induced thermal 
quenching governs disequilibrium CO^CILi chemistry 
depends upon the orbital timescale relative to the mix- 
ing timescale (as a function of K zz ) along the atmo- 
spheric profile. In some cases, vertical quenching along 
the periapse profile may govern the disequilibrium abun- 
dances of CH4 and CO throughout the upper atmo- 
sphere over the entire orbit. For both HAT-P-2b and 
CoRoT-lOb, the effect of orbit-induced thermal quench- 
ing is roughly equivalent to transport-induced quench- 
ing assuming K zz ~ 10 7 cm 2 s _1 . For lower K zz val- 
ues (< 10 7 cm 2 s^ 1 ), thermal quenching may have a 
significant effect, whereas quenching via vertical trans- 
port will determine disequilibrium abundances of CO 
and CH4 at higher K zz values (> 10 7 cm 2 s _1 ). For 
CO« = ^CH4 chemistry, differences in the quenching mech- 
anism may result in large differences in the abundances 
of disequilibrium carbon-bearing species. Refinements 
to atmospheric structure models and improved observa- 
tional estimates of CO and CH4 on highly eccentric ex- 
oplanets may yield important clues to the chemical and 
dynamical behavior of their atmospheres. Further devel- 
opment of the chemical models, including the considera- 
tion of photochemical production and loss rates through- 
out the orbit, may likewise provide improved abundance 
estimates of disequilibrium species throughout the upper 
atmospheres of highly eccentric transiting exoplanets. 
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